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Abstract 


Within the framework of generalized free field theory at nonzero tempera¬ 
ture we address the problem of current conservation. The formalism of thermo 
field dynamics is used to derive a conserved and thermodynamically consis¬ 
tent physical current operator. Consequences for the calculation of photon 
emission rates from a hot plasma are considered briefly. 
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To common wisdom the fundamental interactions of the universe are described by gauge 
theories. Most often they are used perturbatively, i.e., physical quantities of interest are ob¬ 
tained by expansion in terms of non-interacting stable asymptotic states. Such a treatment 
can very successfully account for nearly all data measured in few-body scattering experi¬ 
ments. However, gauge theories also play an important role in physical systems with high 
temperature and a large number of constituents, such as the early universe or relativistic 
heavy-ion collisions. Unfortunately for the theoretical physicist, the perturbative approach 
to gauge theory breaks down at nonzero temperature [|I|. 

Technically this breakdown emerges in the form of severe infrared divergences, which 
one can cure to a certain loop order by the introduction of temperature dependent masses 
for the gauge bosons. The method of “hard thermal loops” has evolved along these lines. 
From a more fundamental point of view however this breakdown is understood through a 
theorem by Narnhofer, Requardt and Thirring ||, which states that (quasi-)particles at 
nonzero temperature must be strictly free of interaction. 

Conversely, to obtain a meaningful solution of a gauge theory at nonzero temperature 
to all orders, one has to deviate from the picture of asymptotically stable excitations. A 
consequent reformulation of a scattering theory in terms of such states with an a-priori 
spectral width, so-called generalized free fields, has been achieved by Landsman |3|. The 
central quantity of this method is the spectral function A, or the imaginary part of the 
retarded Green function S R (E,p) = R eG(E,p) — hr A(E,p) in momentum space. Diagram 
rules, transport theory and other developments of traditional thermal held theory may be 
carried over to this treatment in terms of spectral functions, and one may also show that 
this method is indeed non-perturbative in the sense of not being power-expandable around 
zero coupling constant 0. 

However, particular problems arise with such a formulation, and one of these is addressed 
with the present paper: What is the conserved current of an (abelian) gauge theory in a 
hot system, where the constituents according to the NRT theorem do not have a mass shell 
? We first solve this problem, and then construct the physical current operator for thermal 
held theory. Finally we examine the consequences of current conservation for the production 
rates of gauge bosons (photons). 

To this end we consider a generalized free fermion held in the framework of thermo 
held dynamics (TFD) for quantum helds with continuous mass spectrum |||j. In this 
formulation, the held is described by two different representations ijj x and ij) x , having the 
canonical anticommutation relations 



( 1 ) 


and anticommuting with each other. The statistical doublets of the fermions are, in slight 
modification of ref. 0], introduced as 




where ip x = i^(x). The two-point function therefore is a 2x2 matrix with elements S ab = 
—i /T v k“dR/ ), where (•) = ((1| • | W)) denotes the thermal average (a simple matrix element 
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in the TFD formalism). This propagator equals the two-point function of the Schwinger- 
Keldysh formalism ||, and by construction obeys S ' 11 + S 22 — S 12 — S 21 = 0. Here this 
identity follows from 

((1| (A + ty*7°) = 0 , ((1| (^ x + iqVx) = 0 . (3) 

Since we assume a generalized free field, it is described by a bilinear lagrangian. In contrast 
to ordinary free fields however this lagrangian includes a non-local term, 

£ [T] = T* (r 3 (id„Y -m))^ x - Jd 4 y dH £(ic, y) . (4) 

In this equation, 73 =diag(l,-l) and £ is also a 2 x 2 matrix. Note, that due to the off- 
diagonal matrix elements £ 12 and £ 21 of the self energy function, tildean and non-tildean 
fields are mixed in this lagrangian. 

In £ we now perform a local gauge transformation according to 

-»• (1 - ia x )^x , Vx -> (1 + ia x )ty x , ( 5 ) 

which leads to a variation of the lagrangian 

SC X = §ff ^ 2 , 7 - 37 ^ 2 ; + i jd 4 y iff x (y - xY £(x, y) 

f d 4 y ^*(y~~ +°( d3 ) • ( 6 ) 

After a Fourier transform (we assume a homogeneous equilibrium system), we obtain as the 
conserved current of the generalized free field 

j“(x) = V.TrfV, - i jd‘yd‘z 0- f)t f, A<-(p, % . (7) 

The function A (p,q) occuring here has an expansion in terms of derivatives of the Fourier 
transformed self energy function £(g) 

A-(p.,) = Ae (9)+ ^_^Ai_e (9)+ .... ( 8) 

Current conservation, i.e., vanishing of d^J^(x) on the operatorial level then amounts to the 
requirement 

{p-q)^{p,q) = S(p) -s(g) , (9) 

which is the well-known Ward-Takahashi identity for the matrix valued self energy function. 
It leads us to the identification of 73 7 M — A M (p, q) as the irreducible vertex function of thermal 
field theory. 

I 11 principle, this method may be used to derive the conserved current for any nonlocal 
theory. However, for the special case of thermal field theory with the doublet structure 
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of the lagrangian (j|), it leads to an unphysical current operator: Its thermal groundstate 
expectation value does not depend on the particle density. Note, that this problem is equally 
present in other formulations of thermal held theory 

To construct the conserved physical current which is the source of an external gauge 
held associated with the local phase invariance of the lagrangian ([|), we have to split the 
unphysical operator (|7|) in two parts which are time reversed to each other. The “half” with 
the proper boundary conditions then may be regarded as the physical current operator j^(x) 
of a thermal system with continuous mass spectrum. This splitting is complicated by the 
presence of the mixing terms £ 12 and £ 21 in the lagrangian (3|), 

We find, that the following three requirements are sufficient to determine this operator 
completely: 

1. j fl (x ) is conserved at the operatorial level, dj tl (x)/ dx IJ ' = 0. 


2. j M (x) is a bilinear functional of the fields, and in case of a diagonal self energy function 
with £ 12 = £ 21 = 0 reduces to the if}... i/j component of J M (x), i.e., to “one half” of 
the Noether current. Note, that this requirement also leads to a reduction to the free 
physical current in case of a completely vanishing self energy function, liniy^o = 

^xY^X- 

3. j^(x) has a real expectation value in the thermal ground-state, with a well-defined 
thermodynamical interpretation. 


The operator for the conserved physical current, which follows from these three requirements 

is 


i"(x) = - i 


+ i/ Yy / dr @(x 0 — r) 


^(r,a;)S (r,x;y) 


i7 Yv) s (y\T, X) ^(r, *) 


+ i^° d 4 y drQ(y 0 — r) sign(r/ 0 ) \ *)S 12 (r, x\ y), ip(y) l , 


( 10 ) 


where {-, •} is the anticommutator of the enclosed quantities, and ©(■) is the Heaviside step 
function. 

The construction of the physical current according to a temporal boundary condition, as 
well as the quantization described here, imply the choice of a certain reference frame. Most 
naturally this is the frame where the three-vector current expectation value vanishes. This 
explains the factors in the above operator, which however together with the temporal 
step functions may be generalized to an arbitrary timelike hypersurface. 

The anticommutator term in j^(x) is, due to the fact that we consider generalized free 
fields, a complex number equal to 

7 = ^ JdE (R eG(E,p) - mA(E,p )) A („(E)ImE s (B,p)) . (11) 
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R eG(E,p) and —ttA(E,p) are real and imaginary part of the retarded propagator, and 
Z R (E,p) is the retarded self-energy function. 

The occurrence of this term can be understood when expressing the action of the j lx (x) 
current on the thermal groundstate of the system only through the non-tildean fields. Ac¬ 
cording to eq. (|3|), 


(UK - i7Vw) = {Mi’yl’x 

(Wx (-1^7°) = -«1| 1>ytl> x = (CM tf’x'i’v ~ {Vh> Vv}) • ( 12 ) 


The calculation of the thermal groundstate expectation value of j ,J '(x) current is rather 
tedious and requires a careful consideration of the step functions under the integral. It is 
real and equal to 


(f(x)) = 5»° dE 


Tr 


7° — 


<9R eE(E,p) 
dE 


A(E,P) + aim ^g(f' p) ReG(g.p) 


(13) 


Actually this is identical to the thermodynamical particle density for particles with a nontriv¬ 
ial spectral function 0, i.e., (j 11 (x)) = <5 M ° di1/dp where is the thermodynamical potential 
of the effective held theory and /j is the chemical potential. Hence, all the requirements for 
the physical current operator are fulfilled by our construction. 

Finally we discuss the question where the above derivation might be useful. To this end, 
we assume that our model lagrangian describes only a part of an interacting system, where 
gauge bosons (photons) are coupled to fermions with a nontrivial spectral function. An 
excellent example for such a system might be a strongly interacting plasma of quarks (plus 
other electrically neutral particles). 

In such a system we want to calculate the radiation rate of thermal photons, i.e., the 
“glow” of the plasma [|]]. I 11 principle, this has to be achieved by considering the retarded 
physical current-current correlation function 


nJU® - v ) 


Q(x°-y°) (j»(x)Uy) 


(14) 


If we consider this in momentum space, and fix the coordinate system such that the 3-axis 
is aligned with the photon emission direction, the emission rate of real photons with energy 
(^momentum) k is 0 

R = 2 -y Im (n* + n”) (15) 

= 2 Im (-nyj only if k * n,„i/'i = 0 . (16) 

In principle, this calculation would be possible given our conserved current (0) - but the 
result would involve the irreducible vertex function to a given order, and hence might be 
difficult to obtain. 
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One may therefore ask, in which cases it is sufficient to calculate the retarded one-loop 
fermion polarization tensor for our nonlocal effective held theory, 


= -i 


• d a p 

Wr 


Tr [a In S R (p + k)g 7„ S'* (po, p)] 


(17) 


which differs from the correlation function in two aspects: It does not have the proper vertex 
factor (which is of higher loop order), and it also misses the (^-corrections of (|7|). 

Clearly this is the case if and only if the retarded self energy function E R (E,p) is a func¬ 
tion of the energy E alone (and independent of the three-momentum p of the particle). The 
other matrix elements of E are obtained from the retarded self energy by a transformation 
depending only on this energy parameter [|], and hence the vertex correction according to 
eq. (j8j) is also proportional to 

Independence of E^ from the momentum might be good approximation especially for 
fermions coupled to massless bosons 


10 , and it follows that 


dE R (E,p) 

dp 


o^ng ]R (fc) = ng(fc), i,j = 1,2,3; 


(18) 


whereas ^ Tiff (k). Hence, for systems approximated by a momentum-independent 

retarded self energy function, one may calculate the photon radiation rate from the effective 
one-loop polarization tensor by use of eq. (p~5|) , but not by use of eq. (|T6| ) j8j. 

With the present paper we reached two goals. First of all, we have established a clean 
path towards current-conserving (and thus gauge invariant) calculations in generalized free 
held theory. As we have shown, this is especially relevant for the determination of emis¬ 
sion rates of gauge bosons (photons) from a hot plasma. Secondly, we have constructed a 
conserved physical current operator in the framework of thermo held dynamics, which is 
thermodynamically correct. To our knowledge, this has not been achieved in other formu¬ 
lations of thermal held theory, we are currently exploring this connection in more detail. 
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